Controlling the generation of high frequency electromagnetic pulses with relativistic 

flying mirrors using an inhomogeneous plasma 
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A method for the controlled generation of intense high frequency electromagnetic fields by a break- 
ing Langmuir wave (relativistic flying mirrors) in a gradually inhomogeneous plasma is proposed. 
The wave breaking threshold depends on the local plasma density gradient. Compression, chirp- 
ing and frequency multiplication of an electromagnetic wave reflected from relativistic mirrors is 
demonstrated using Particle-In-Cell simulations. Adjusting the shape of the density profile enables 
control of the reflected light properties. 
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High frequency coherent radiation sources are of great 
interest for various applications and future understand- 
ing in fundamental science [H, 0. In material science, 
the study of extremely short time reactions requires high 
temporal and spatial resolution imaging which could be 
achieved using extremely short duration high frequency 
electromagnetic pulses. High-power X-ray radiation can 
be used to analyse nano-scale objects such as complex 
molecules offering low noise and sufficient energy be- 
fore their destruction due to Coulomb repulsion [3[. In 
laser-matter interaction, the possibility to develop high 
frequency extremely-intense laser systems would lead to 
novel quantum electrodynamics effects Q. 

An approach to generate high frequency radiation 
based on the concept of the flying mirror considers a 
plasma shell travelling close to the speed of light as 
a relativistic mirror. Reflected light undergoes double 
Doppler frequency up-shift, compression, intensification 
and focusing due to relativistic effects. Various schemes 
were described 043 an( i experimentally demonstrated 
Q as a proof of the feasibility of this concept. 

Controlling the generation time and position of the 
breaking wave is an essential aspect for experiments to 
synchronize the mirrors with an incident laser pulse in 
an accurate, flexible and facilitated way. It also allows 
for tuning the reflected pulse properties. 

In the laser wakefield acceleration regime wave 
breaking offers electrons the initial injection required to 
enter the acceleration phase of the wakefield [l0(. For 
this purpose, tailored inhomogeneous plasmas have been 
considered as a way to control the wave breaking. In a 
downward density profile, the plasma frequency depends 
on the coordinates. In the wave, electrons progressively 
oscillate out of phase resulting in the breaking of the 
wave [ll|. This mechanism was examined in Refs. 0- 



* Also at Prokhorov Institute of General Physics, Russian 
Academy of Sciences, Moscow 119991, Russia 



[l4[ in the case of an inhomogeneous plasma with scale 
length larger than the plasma wavelength. Step-like 
density transitions are also useful for controlled electron 
injection [IlGjl. They can be produced experimentally 
using a transverse laser pulse [Hf or using an obstacle 
(razor blade) inserted into the ultrasonic gas flow [l7j . 
Recent experiments in Rcfs. [HI use a density deple- 
tion channel generated in a gas-jet by a transverse laser 
pulse. 

Under the realistic conditions of experiments, the 
plasma is always inhomogeneous. In the present paper, 
we examine plasma inhomogencity effects on the flying 
mirrors and propose using plasma inhomogencity as a 
new way to control the generation of localized flying 
mirrors. We use ID Particle-In-Cell (PIC) simulation 
for detailed analysis of the formation of the mirror and 
the reflection of the light. The possibility to tune the 
properties of the reflected light by adjusting the plasma 
parameters is demonstrated. 

An electromagnetic pulse travelling in an undcrdense 
plasma (i.e. the local density n is below the critical den- 
sity n c — m e ui 2 / Aire 2 where ui is the laser frequency, m e 
is the mass of an electron and e its charge) generates a 
driven electron plasma wave. In the linear regime, the 
wave has a sinusoidal shape and electrons osc illate at the 
non-relativistic plasma frequency lo p = \J Ann e e 2 /m e 
where n e is the local electron density. The plasma wave 
propagates with phase velocity close to the group veloc- 
ity of the laser driver pulse v p h — v g Q . 

For a sufficiently intense driver pulse, one enters 
the nonlinear regime. The plasma wave breaks when 
Vph < v e with v e being the electron velocity. Nonlinear 
effects cause the steepening of the wave. The electron 
density appears strongly modulated with high de nsity 
shells formed at the positions of maximum velocity [19| . 

In an undcrdense plasma, the Lorentz factor calcu- 
lated for v e , 7 e = 1/ \/\ — (v e /c) 2 , can be approximated 
by 7 e = 1 + a 2 /2 assuming a stationary wave solution, 
where a = eE/mcuj is the normalized laser amplitude 
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FIG. 1. Tailored inhomogeneous density profile in the 
framework of the flying mirror. 



wave. It yields k p 



«p,0 



d x L)pt. The phase velocity, 



equal to v p h = uj p /k p , therefore is given by 
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where k p $ is the initial wave number. The plasma fre- 
quency uj p depends on the local value of the electron 
density n e (x) which yields d x to p = yj ire 2 /m e d x n e / ^JrQ. 
In a gradual inhomogeneous density profile I^Wpl -C 1 
we have 



/ 7re 2 d x n e 

Vph - Vphfi 1 + \ 1 
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and E is the laser electric field. The Lorcntz factor 
for the phase velocity, ^ p h = 1/ yl — (fp'i/c) 2 , is re- 
lated to the plasma frequency as j p h = (oj/u p )\/1 + a 2 , 
where the dependence of the laser group velocity v g on 
the laser amplitude has been taken into account. Writ- 
ing the breaking condition 7 e = 7^ [T9j in the form 
a 2 / 2 + 1 = uj/uj p \/l + a 2 , we obtain [2(j the breaking 
threshold 



_ f ^2 (w/wp - 1) for 
(2w/ Wp ) 2 / 3 for 



a < 1 
a > 1 



(1) 



Strong density modulations formed in the wake of an 
intense driver pulse constitute relativistic flying mirrors. 
A counter-propagating electromagnetic pulse, referred 
to as the source pulse, is partially reflected back with 
frequency up-shifting and compression. For normal in- 
cidence reflection, the reflected pulse frequency lo t is 
given by the relationship 



l + Pph(t) 



(2) 



for P p h(t) = v p h{t)/c where changes in (3 p h{t) are much 



slower than 



In the case of a relativistic mirror 



ph 



1, Eq. ([5]) can be approximated by u> r ~ 47, 



ph u 



The reflected electromagnetic field E r is amplified by 
the same factor so that E r = E s (l + f3 p h(t))/(l — (3 p h(t)). 
In the general case, the dependence on time of (3 p h re- 
sults in the chirp of the reflected light [2l[ . 

The reflected pulse energy is determined by the mir- 
ror reflectivity. For a small amplitude plasma wave, the 
reflectivity is exponentially small. Close to the wave 
breaking, the electron density becomes singular provid- 
ing high reflectivity with the number of reflected pho- 
tons being proportional to 1/ j ph [22j ■ 

In an inhomogeneous plasma, the plasma wave be- 
comes a wave with a continuous spectrum [ill . fl2T | . Its 
wave number, k pi depends on time t and position x ac- 
cording to the equation 



d t k 
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This equation is a consequence of the relationships uj = 
—dfO, k = d x 9 and the cross differentiation property 
d tx 8 = d xt where 6(x,t) is the eikonal of the plasma 



where v p h,o is the phase velocity in a homogeneous 
plasma. In a downward density profile d x n e < 0, the 
phase velocity progressively decreases with time. The 
corresponding deceleration depends on the magnitude 
of the density gradient |9 x n e |. The breaking condition 
Iph = 7e is inevitably reached after some wave periods. 
The breaking time % r [l3| is equal to 



2n K 
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In order to study the flying mirror properties in inho- 
mogeneous plasmas we perform numerical simulations 
using a ID 2/2 parallel PIC code. 

Fig. [T] illustrates the assumed tailored density pro- 
file composed of two homogeneous plasma slabs of den- 
sity m = 0.025n c (~ 0.2 75 x 10 20 cm" 3 x (I^m/A) 2 ) 
and 712 = 0.0175n c corresponding to plasma wave- 
lengths of Ap ei i ~ 40A and A pe! 2 — 47. 5A and re- 
spective lengths L\ — 200A and L2 = 350A where 
A = cuj^ 1 . They are separated by a decreasing inhomo- 
geneous region having a cosine square density profile, 
n e (x) = (ni — n-i) cos 2 (irx / L) + TI2, which provides con- 
tinuous density variation from n\ to 712 with the length 
L = 150 A (~ 24 /im for A = 1 /mi) which approximately 
corresponds to three plasma wavelengths. To limit the 
effect of the vacuum heating [23j , a pre-plasma of length 
Lq = 200coj~ 1 gradually rises from zero to n\ at the left- 
side of vacuum-plasma interface. 

The driver pulse is Gaussian linearly polarized along 
the y direction with wavelength A^ = A, longitudinal 
full width at half maximum (FWHM) 1<j = I0A and its 
amplitude is equal to a^ = 1.5 (/ ~ 3.1 x 10 18 Wcm~ 2 x 
(1 /zm/A) 2 ). Densities in the homogeneous regions give 
for the breaking condition a& r = 3.2, from Eq. [I] so that 
the generated wake wave is initially below the breaking 
level ad < ab r - The source pulse is rectangular and is 
polarized in the direction perpendicular to the driver 
polarization. Its length is equal to l s = 80A, X s — 2A 
and the amplitude a s = 2 x 10~ 4 is relatively low to not 
disturb the driven plasma wave. 

In Fig. [2ji, the simulation results are presented. The 
driver pulse propagating in the positive x direction ini- 
tially generates a weakly nonlinear plasma wave. The 



3 



0.25 



,0.00 



-0.25 



Inhomogeneous region 



> 1 


i 
i 


a) 


— iHr^ 


L J 1 






1 Z | 




\ Reflected pulse | 


Driver pulse 



900 



1000 1100 1200 




b) 



1050 



- ut 
oXlO 1500 1550 1600 1650 1700 1750 
K) 2.5F 




o 1 — 1 



10' 




/ Px,c 


e) 




w 

1 1 



1530 1550 1570 

LOt 



970 990 1010 1030 



FIG. 2. Results of the PIC simulation: a) Generation of the 
density spikes in the plasma wave (n e ) and reflection of the 
source pulse (E z ). b) Evolution of the plasma wave phase ve- 
locity P v h in the inhomogeneous plasma region. Solid curves 
correspond to the theoretically calculated j3 p h, markers give 
P P h found in the simulation and the dotted line indicates the 
position of the minimal j3 p h- c) Time- frequency spectrum 
|i? z (u;,£)| of the electric field E z (lu) from the light reflected 
by the mirrors, d) Analysis of a parabolic chirp. The line 
corresponds to the analytical evolution of the frequency, e) 
Acceleration of electrons by the breaking plasma wave. A 
low pass filter was used to extract the electric field of the 
reflected light from the main source pulse. 



wave only breaks in the inhomogeneous region. Density 
spikes develop in a finite time and constitute temporary 
existing mirrors. Since we know the breaking time (|7J|), 
this process ensures spatial and temporal control of the 
breaking wave. 

The source pulse is longer than the plasma wave- 
length, the reflection on the periodic train of density 
spikes generates a complex signal composed of several 
wave packets. 

The time dependence of the wake wave phase veloc- 



ity versus position in the inhomogeneity is shown in Fig. 
[2Jd for different times. The numerical results (shown by 
markers) appear in good agreement with the theoret- 
ically decreasing evolution of the phase velocity (solid 
lines) given by Eq (0}. The breaking condition is at 
first satisfied at t = lSOu;^ 1 after the driver pulse has 
entered the target. The position of the minimal phase 
velocity progressively tends to the position of the max- 
imal density gradient. 

The time-frequency spectrum of the reflected radia- 
tion is presented in Fig. where the electric field E z is 
plotted in the plane (uj, t). The time-frequency analysis 
of E z (t) is performed using the Gabor transformation 
defined as 12411 
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where a = 41og(2)c 2 /^ so that l w is the full length at 
half maximum of the Gaussian window. 

A remarkable property of the high frequency radiation 
is the generation of a train of short electromagnetic wave 
packets. Space and time evolution of the mirror velocity, 
v p h, leads to the generation of parabolic chirps as shown 
in Fig. Hi. 

In order to calculate the frequency chirp we consider 
a plasma frequency, in the inhomogeneous region, of the 
form 



uj p (x) 



Auj 



/ x 

— tanh — 
2 \L 



LUl + U)2 



(S) 



with Aw = cji — L02- From Eq. @, the phase velocity 
can be locally approximated by 



Vph Vph.rnin "t - &X 
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around the position of minimal phase velocity v p h_ rnin 
with 8 = Aujt/2L 3 k p . Using Eq. we obtain for the 
frequency multiplication the following quadratic equa- 
tion 



Ld r — CJ yVl 



i>ii 

.nit!) 



(10) 



where Lu r ,min is the frequency multiplication after re- 
flection on a mirror travelling at Vph.min • 

Eq. (nnj) 

characterises the parabolic frequency modulation of the 
reflected pulses. 

When electrons are injected into the wakeficld, the 
mirrors decelerate due to momentum conservation 
which causes the wave packets to be down chirped. 
This phenomenon happens at the reflection from the 
fourth mirror formed behind the driver pulse as shown 
in Fig. [2};. Due to the decreasing mirror velocity v p h, 
density maxima gradually increase offering higher re- 
flectivity. The life spans of the mirrors are gradually 
prolonged. The average reflected frequency diminishes 
and the spectrum broadens according to Eq. (|2|). 

The properties of the reflected pulse can be modi- 
fied by changing the parameters of the inhomogeneous 
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FIG. 3. a) Number of reflected photons normalized by the 
number of incident ones N p h, r /N p h,i- b) Average reflected 
frequency for density variation An e changing from 0.005n c 
to 0.0175n c as a function of the plasma wave cycles at which 
the reflections have occurred. 
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FIG. 4. a), b), c) Normalized density n e and electric field 
E z at time t = 1100w _1 for L = 70, 50 and 30X. d) Initial 
normalized density pit n e . e) Spectra |Z?~(aj)| of the reflected 
wave packets for the different pit lengths. 



region: density profile, length L and density variation 
An e = m — n<i. 

In the spectral domain, average reflected frequency 
can be defined as 



\E(u>)\u>du> 



(11) 



\E(u))\du 



and the number of photons as 
50 E(uj) 2 



N ph = 



8ir 2 huj 



du3 
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For density variation An e changing from 0.005n c to 
0.0175n c , the evolution of the reflected number of pho- 
tons divided by the numbers of incident ones is shown 
in Fig. [3^ as a function of the plasma wave cycle when 
reflection occurred. The average frequency, Q r , is plot- 
ted in Fig. [3Jd. The length of the inhomogeneous re- 
gion is equal to L = 200A. We interpret the increase of 
the reflected photon number as the onset of the wave 



breaking. As predicted by Eq. ©, the larger is the 
density variation An e , the shorter the breaking time. 
For An e = 0.005n c , the wave breaks at 1380J" 1 af- 
ter the driver pulse has entered the inhomogeneity re- 
gion, three wave periods behind the driver pulse. For 
An e = 0.020n c , wave breaks during the first plasma 
wave cycle. For a given density, the frequency decreases 
while the number of reflected photons grows. Progres- 
sive decreasing of the phase velocity gives higher reflec- 
tivity but lower frequency up-shift (Eq. [2]). 

Now, we consider a localized density depression, as re- 
alized experimentally in [l8| , with a scale length shorter 
than the plasma wavelength L < A p .i- The density de- 
pression profile is composed of a first decreasing plasma 
region of length L followed by a symmetric increas- 
ing one. As illustrated in Fig. @Ji, the density varies 
from n e = 0.005n c to a minimal value n e = 0.00 457i c 
(A p ~ 90A) for lengths L = 30, 50 and 70A. The nor- 
malized amplitude of the driver pulse is equal to a c i = 3 
whereas other parameters remain the same as previously 
described. 

In the short-scale density perturbation, wave break- 
ing leads to the generation of very short-lived mirrors. 
Reflected pulses have a very short duration (~ 3 fs for 
A = 1 fxm), composed of a few cycles as displayed on Fig. 
[4ji, b and c. Decreasing the depression length reduces 
the reflection time and generates shorter wave packets. 
For the density depression lengths L = 30A (Fig. |4j:) 
and 50A (Fig. 0Jd), the density gradient is sufficiently 
abrupt to result in the acceleration of an electron bunch 
seen in the density profile as a narrow density spike. 
This is of interest for the electron injection scheme fl2T - 
[Hj]. In this case, a characteristic down chirping of 
the reflected light can be observed. The spectra of the 
reflected signals in Fig. 2J> show that the frequency was 
multiplied by a factor of 16 for length L = 70A. The 
shortening of the density depression length results in the 
generation of a stronger negative chirp. 

In the downward density region, the plasma wave pro- 
gressively travels with lower velocity than the reflected 
light. In the rising density part of the density depres- 
sion, the plasma wave accelerates so that the distance 
to the reflected pulse can be reduced. In the simula- 
tions, the reflected pulse is not affected by the density 
modulation in the rising inhomogeneous density region. 

In conclusion, a tailored downward plasma density 
offers time and space control for the generation of rela- 
tivistic mirrors in a breaking laser-driven plasma wave. 
The reflected light properties including frequency mul- 
tiplication, compression, chirping and duration can be 
modified by adjusting the density inhomogeneity pa- 
rameters. This discovery is useful for experimental stud- 
ies on the flying mirror concept, including compression 
of controllably chirped attosccond pulses with aperiodic 
soft and hard X-ray mirrors [2o| . It makes a new step 
in the realization of a stable, coherent, compact and 
tunable high frequency light source required for various 
applications. The analysis of the reflection of the laser 
source on a density profile can also be contemplated as 
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a plasma diagnostic tool. 
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